Orbital-selective Mott transitions in a doped two-band Hubbard model 

with crystal field splitting 
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We investigate the effects of crystal field splitting in a doped two-band Hubbard model with 
different bandwidths within dynamical mean-field theory (DMFT), using a quantum Monte Carlo 
impurity solver. In addition to an orbital-selective Mott phase (OSMP) of the narrow band, which 
is adiabatically connected with the well-studied OSMP in the half-filled case without crystal field 
splitting, we find, for sufficiently strong interaction and a suitable crystal field, also an OSMP of 
the wide band. We establish the phase diagram (in the absence of magnetic or orbital order) at 
moderate doping as a function of interaction strength and crystal field splitting and show that also 
the wide-band OSMP is associated with non-Fermi-liquid behavior in the case of Ising type Hund 
rule couplings. Our numerical results are supplemented by analytical strong-coupling studies of spin 
order and spectral functions at integer filling. 

PACS numbers; 71.10.Fd, 71.10.Hf, 71.30.+h, 71.27.+a 



I. INTRODUCTION 

As a fundamentally nonperturbative phenomenon be- 
yond simple electronic band pictures, the Mott-Hubbard 
metal-insulator transition has been a subject of great 
interest in solid-state physics for decadesj^ In recent 
years, the Mott-Hubbard transition drew significant ad- 
ditional attention in the context of ultracold atoms. ^'^ 
In the traditional scenario, merging the suggestions by 
Hubbard^ and Brinkmann and Ricej^ the electrons in 
a half-filled valence band acquire more and more effec- 
tive mass with increasing interactions until they local- 
ize simultaneously and form a paramagnetic insulating 
state; this picture has found support in numerous calcula- 
tions within dynamical mean-field theory®'^ (DMFT) for 
the Hubbard model^i^i^ and its extensions to multipleiS 
equivalent (e.g., to three t2g) orbitals. However, more re- 
cent experiments on ruthenate a^^'^^ have indicated that 
interesting complications can arise in the presence of 
multiple inequivalent valence orbitals with different ef- 
fective bandwidths: then, as first illustrated in a sim- 
plistic model without inter-orbital coupling, increas- 
ing interactions U > Uci first localize the electrons in 
the narrow orbital(s) while the wide-band electrons ini- 
tially remain itinerant and become insulating only at 
U > Uc2 > Uci- Subsequently, such orbital-selective 
Mott transitions (OSMTs) have been discussed for vari- 
ous classes of materials;-'^'^ the idea of partial localization 
has also been extended to "momentum-selective Mott 
transitions" in the context of high-Tc materials J^ii^ 

While it was immediately clear that a realistic de- 
scription would require at least three bands for the 
ruthenates^^ and that the mechanisms leading to differ- 
ent band widths would generically also affect the band 
center (i.e., induce a crystal field), most studieaii"— so 
far have addressed a minimal 2-band Hubbard model 



(1) 



in which the two orbitals m S {1,2} differ only by their 
hopping amplitude (between nearest-neighbor sites 
i, j) and share the same local intra-orbital interaction U 
and site potential; also the third term, containing both 
an inter-orbital Coulomb repulsion, parametrized by U' 
(with < U' < U), and an Ising- type Hund rule cou- 
pling with amplitude Jz > is symmetric in the orbital 
index. The DMFT studies of this model have almost 
exclusively assumed a semi-elliptic "Bethe" density of 
states and the absence of any magnetic order and mostly 
focussed on half fiUing {n = 2) and (to a lesser degree) on 
a hopping ratio of ^1/^2 = 2. The expectation that this 
"standard model" for OSMTs captures the essence of the 
phenomenon, i.e., resolves two distinct orbital-selective 
transitions at half filling and low temperatures, could in- 
deed be verified, ^■^ after some initial confusiouj^l In ac- 
cordance with previous literature) ^^'^^ we will refer to 
this model as the "J^-model". The general doped case 
was investigated in detail in Ref. [l^ with the use of quan- 
tum Monte Carlo (QMC) methods and previously also by 
exact diagonalization^ and slave boson methodsi^ 

In the following, we will explore the more general 
Hamiltonian 
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Here, the third term, proportional to J±, describes spin 
fiips and pair hopping processes arising from the general 
Hund rule coupling. The second term, proportional to 
A, shifts the relative positions of the atomic energy lev- 
els of the two orbitals and, hence, describes crystal field 
splitting. First results regarding the impact of the lat- 
ter term in the three-band extension of Ha have been 
obtained within a slave boson approximation.'^" Other 
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studies have included crystal field terms as in ([2]), but for 
orbitals with identical band widths<^"— Our goal in this 
paper is to explore the physics of TJa (with 7^ t2 and 
A 0) using QMC simulations within DMFT, identify 
new phases, and discuss spectral properties. An impor- 
tant special case of H/\ occurs for J± = Jz = J, which 
we refer to as the " J-model" . 

The J-model with ti = t2 and A = is rotationally 
invariant in the sense of Castellani et al<^ We note, how- 
ever, that the experimental systems allegedly described 
by H^, including the perovskite ruthenates, have less 
than cubic symmetry, so that there is no reason for as- 
suming a rotationally symmetric screened Coulomb in- 
teraction in ([2|) , if the asymmetric hopping and the spin- 
orbit interaction are taken into account. Crystal field 
splitting of d-orbitals also reflects a broken rotational 
symmetry. For this physical reason, we will below set 
Jj_ = in quantum Monte Carlo (QMC) calculations. In 
analytical arguments we will assume < Jl < Jz , which 
allows for antiferromagnetism at strong coupling. The 
choice Jl = in numerical simulations has the impor- 
tant additional advantage of avoiding a sign problem. 

Note that, irrespective of the value of J±, our Hamil- 
tonian is particle-hole symmetric only for A = 0; for 
A ^ it is mapped to i?_A under a particle-hole trans- 
formation. Accordingly, we can cover the entire density- 
and A-range by assuming n > 2 and A € K; results for 
n < 2 then follow from particle-hole symmetry. 

This paper is built up as follows. First, in Sec. |IT1 we 
consider the analytical properties of the Jz- and the J- 
models. In particular, we calculate the strong-coupling 
Hamiltonian and the local spectral functions at half fill- 
ing; we also comment on the strong-coupling Hamiltonian 
at quarter filling. Then, in Sec. lIIIl we discuss the results 
of our QMC simulations for the Jz-model, in particular 
for orbital occupations, spectral functions and for the 
phase diagram. We also comment on non-Fermi-liquid 
properties on the basis of imaginary-frequency self-energy 
data. Generally, for the QMC calculations, we will con- 
centrate on the question of metallicity as effect of cor- 
relations; hence, we restrict ourselves to paramagnetic 
phases. We close (in Sec. ITV)) with a summary and an 
outlook. Technical details regarding the strong-coupling 
expansion of the spectral functions are deferred to an 
appendix. 



at A = and ti = t2, it is also rotationally symmetric 
(i.e., S'0(2)-symmetric) in the band index. In addition, 
the model with A = is particle-hole symmetric at half 
filling. Crystal field splitting breaks the rotational sym- 
metry for J± = Jz and the particle-hole symmetry at 
half filling. 



A. Strong-coupling Hamiltonian at half filling 

For the half-filled model H/^ with J± — Jz and 
A = 0, an effective strong-coupling Hamiltonian was de- 
rived by Ferrero et al3S Since most of their arguments 
are also valid for the model without particle-hole sym- 
metry (A ^ 0), and calculations for the model with 
J± < Jz are very similar, we only sketch the deriva- 
tion. The effective strong-coupling Hamiltonian is ob- 
tained from standard Harris-Lange degenerate perturba- 
tion theory, which is based on a canonical transforma- 
tion from the Hubbard electrons Cima to new particles 
CimtT, whose hopping leaves the associated total number 
of double occupancies invariant (see also App. For 
these new particles, we then define annihilation opera- 
tors of double occupancies di^ = Ci2a-Ciia- (with a =tii) 
and dio = ^(ci2;Cii^ -|- Ci2^Cii;), and a (5 = l)-spin, 
constructed from these dio-- and dio-operators: 



Si 



dt^d.„ 



dl d 



dt, d 



The relevance of the di^-- and diQ-operators is that the 
three atomic states dj^lO), dj^lO) and dig|0) are lowest 
in energy (and, hence, span a triplet) for J± = Jz = J, 
while the two states d;^|0) and d!, |0) are lowest in energy 
if < Jl < Jz. With these definitions, the effective 
S = I spin Hamiltonian for Jj^ Jz = J is given by: 



= Z^'^Hcis(Si • Sj -Jiinj) , Jjj^i, 



U + J 



where the number operator is defined as rii = '^icr'^i<T^" 
djpdjQ. The effective Hamiltonian for < Jj^ < Jz reads: 



II. ANALYTICAL RESULTS FOR THE J,- 
J-MODELS 



AND 



In this section we present analytical results for the Jz- 
and J-models with general A ^ for integer fillings at 
strong coupling {U -^00). Results of interest include the 
effective strong-coupling Hamiltonians at half- and quar- 
ter filling, which provide information about the symme- 
tries and low-temperature phases of the models, and the 
spectral functions. We note that the model Hamiltonian 
-ZJa is S'[/(2)-symmetric in the spin sector for J± = Jz] 



iJj =^ Ji,(S'i3 5j3-ninj) , Jj^ = 

(ij> 



(^1)^ + it2r 
U + Jz 



where now ni = J2a 'Ha'^ia- N^ote that the value of J± be- 
comes irrelevant if Jz > J_\_. These results are generally 
valid in any spatial dimension. 

It is interesting to note that the results stated above 
for the effective Hamiltonians at strong coupling are en- 
tirely independent of the crystal field splitting parameter 
A, although they are valid for general A <[/ -I- Jz if 
J± < Jz- The explanation is, technically, that the sum 
of the excitation energies for hopping processes, starting 
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from the subspace spanned by dj|.|0), dl^\0) and possi- 
bly djglO), is independent of A. This, in turn, occurs 
since the crystal field splitting term in the Hamiltonian 
commutes with the hopping. 

In the limit of high spatial dimensions, the spin Hamil- 
tonians, thus derived, are solved by mean field theory. In 
particular, the "antiferromagnetic" critical temperatures 
for these models can easily be calculated for a bipartite 
lattice in the limit of infinite coordination number Z: 



for J± = Jz 
for Jj^ < Jz 



J 



These results are interesting, because they show that 
is dominated in both models by the largest hopping am- 



plitude, which is that of the broad band (^2 = t2\/Z). 
We conclude that the broad band primarily determines 
the energy scale at strong coupling; from previous workSS 
we know that the energy scales of antiferromagnetism at 
weak coupling are primarily determined by the narrow 
band {tl = t\\fz\ Note that, since the effective strong- 
coupling Hamiltonians are A-independent (within their 
range of validity), the same holds for the critical temper- 
atures. 



B. The spectral functions at half filling 

The calculation of the spectral functions on the Bethe 
lattice (which is of interest here) at strong coupling, 
T = 0, and half filling proceeds along the lines of Kali- 
nowski and Gebhard.— These authors calculated the 
spectral function for a sm^fe-band Hubbard model in 
strong-coupling perturbation theory. The analogous cal- 
culations for the two-band model iJA, considered here, 
are technically considerably more complicated and are, 
hence, summarized in App. [X] The results are, in terms 
of the noninteracting density of states z^J^(a;): 



,LHB 



,LHB 



1 



U+JzT^ 



U+Jt^ 



{J± < Jz) 
iJ± = Jz] 



Here the upper sign refers to the (to = 1)- and the lower 
sign to the (to = 2)-orbital. It is interesting to note that 
the spectral function for J± ~ Jz = J is broader than 
that for J± < by a factor of ^3/2, which reflects 
the larger number of possible hopping processes in the 
former case as compared to the latter. Clearly, if Jz — 
J_L is positive but small, there will be a crossover from 
the spectrum for Jz > Ji_ to that for Jz = J± at finite 
temperatures (T > 0) or finite Hubbard interaction ([/ < 
oo). The relevant temperature and [/""'^-scales are (Jz — 
Ji^/k-Q and (Jz — J±)l{ti^2f , respectively. 

We note that the effect of crystal field splitting at 
strong coupling is simply to shift the (to = l)-band en- 
ergetically to the right and the (to = 2)-band to the 



left. We will see below from the results of the QMC- 
simulations, that the effect of crystal field splitting at 
finite interaction strength is not quite so obvious. 



C. Strong-coupling Hamiltonian at quarter filling 

The determination of the strong-coupling Hamiltonian 
at quarter or three-quarter filling is extremely simple, al- 
though the result is nontrivial. E.g., in the atomic limit 
(ti = ^2 = 0), the minimal energy at quarter filling is 
obtained for an ensemble of singly-occupied sites (if |Jz| 
and |J_l| are both smaller than U and U'). The only 
relevant term left in the Hamiltonian is, therefore, the 

crystal field splitting \^^icri''^iia ""-120-)' implying that 
the ground state occurs for an ensemble of single occu- 
pancies in band 2 if A > or, alternatively, single oc- 
cupancies in band 1 if A < 0. The spins of these sin- 
gle occupancies are not fixed yet. Since all particles at 
quarter filling occupy one single band, we conclude that, 
if a small hopping of the particles is now switched on, 
the two-band Hamiltonian at quarter filling reduces 
to a half-filled single-band Hubbard model. At strong 
coupling, the half-filled single-band Hubbard model re- 
duces to an antiferromagnetic Heisenberg model. Hence 
we conclude that, at low temperatures, the quarter-filled 
two-band model JTa describes an orbitally ferromagnetic 
Heisenberg spin-antiferromagnet. 

We note that TJa at three-quarter filling is mapped by 
a particle-hole transformation to 7J_a at quarter filling; 
the low-temperature physics is, again, that of an orbitally 
ferromagnetic Heisenberg spin-antiferromagnet. 

These strong-coupling arguments assume that the 
bandwidths Wi and W2 of the two-band model are small 
compared to all other parameters, so that, in particu- 
lar, |A| > Wi^2- If this condition is not fulfilled, the 
results may change completely. For instance, in the spe- 
cial case A = the two-band model at quarter filling 
changes its behavior both in the spin and in the orbital 
degrees of freedom and reduces to a spin-ferromagnetic 
Heisenberg orbital-antiferromagnetj^ The numerical re- 
sults show that the system, in the ground state, may or 
may not have orbital long-range order, depending on the 
model parameters. 



III. QMC RESULTS FOR THE J^-MODEL 

In this section, we present QMC results for the Jz- 
model ([T]) with Jz = J7/4 and U' = U/2. We will assume 
half-elliptic "Bethe" densities of states for both bands, 
with a full bandwidth Wn = 2 for the "narrow" band 
and Ww = 4 for the "wide" band, respectively; this cor- 
responds to the best studied case for A = 0. Also in 
line with earlier work, we will restrict ourselves to the 
paramagnetic case, i.e., exclude antiferromagnetic and 
orbital order; the results should be relevant at interme- 
diate temperatures or for frustrated systems. In the 
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A (crystal field splitting) 

FIG. 1: (Color online) Orbital occupation numbers Un, 
as a function of crystal field splitting A: (a) for total filling 
n = 2.1 and (b) for n = 2.2. Results for the wide (narrow) 
band are indicated by thick lines and symbols (thin lines). 
Mott plateaus in n-n, appear for U > 2 (see magnified 
insets for riw). 

following, we present QMC results obtained at a temper- 
ature T — 1/40 with a discretization parameter in the 
Trotter decomposition'^'^ of Ar = 0.4; these data were 
checked by additional calculations at different values of 
Ar. 

In the following, we will first discuss orbital-specific 
occupation numbers and spectral functions, then con- 
struct the phase diagram and, finally, address possible 
non-Fermi-liquid properties on the basis of self-energy 
data, all as a function of crystal field splitting. 

A. Orbital occupation numbers 

As the crystal field splitting A acts like a magnetic field 
in the orbital sector, i.e., shifts all energy levels of the 
wide band downwards for A > and all energy levels of 
the narrow band upwards, one expects that, generically, 
an increase of A will increase the filling in the wide band 
and decrease the filling in the narrow band. This applies 
both at constant chemical potential ^ and at constant to- 
tal filling n = nn-l-n„, unless one or both of the bands are 
incompressible, i.e., in a Mott or band insulating state. 



Only at integer filling, both bands can be incompress- 
ible at the same time: indeed, the model is known to 
be fully insulating at n = 2, [/ > 2.8, and A = Oj^ 
due to the first-order character of the associated tran- 
sition, the insulating state must be stable also at small 
A (up to the order of the gap in the particle-hole sym- 
metric case). At very large crystal field splitting (|A| > 
max{Wv!,Wn,U, kB T}), one expects that both orbitals 
become band-insulating, with all electrons in the lower 
band (wide or narrow, depending on the sign of A) . 

Away from half filling (here we concentrate on n > 2), 
at least one of the bands must have a noninteger occupa- 
tion and, therefore, remain metallic. This case appears 
much more interesting, so we will focus on it in the fol- 
lowing. In Fig.[IJa), we present data for n — 2.1, i.e., for 
relatively small electron doping. The case of larger dop- 
ing will be discussed later [cf. Fig. [IJb)] while results for 
hole doping [n < 2) follow from symmetry. As expected, 
the occupation n„ of the wide band (symbols and solid 
lines) increases monotonically with A at all interactions 
2 < U < 3.6. Accordingly, the occupation = n — n„ 
decreases monotonically. 

Plateaus in the orbital occupations are observed for 
large absolute values of the crystal field splitting at all 
interactions beyond thresholds that increase with U, e.g., 
for A > 4.0 and for A < —3.8 aiU = 3.6 (squares), which 
correspond to a totally filled wide and narrow band, re- 
spectively. In between these (orbital selective) band in- 
sulating phases, the density curve ^^(A) is smooth with 
strictly positive slope at [/ < 2.0 (diamonds), indicating 
a purely metallic phase. 

At U = 2.4 (triangles), an additional plateau appears 
at < A < 0.8. On this plateau, the occupation of the 
narrow band is integer (71,1 — 1.0) and that of the wide 
band fractional (tIw — l-l)- Since A acts as an orbital- 
dependent chemical potential (in addition to ^), we can 
interpret the pinning of the narrow-band occupation at 
half filling as a arising from an incompressibility of this 
subsystem; thus, the system appears as a narrow-band 
orbital-selective Mott insulator. In contrast, the "pin- 
ning" of at the value 1.1 arises from our constraint 
?i = rtn + = 2.1; it would not show up in plots at con- 
stant /i. The narrow-band Mott plateau broadens as the 
interaction is increased and extends to slightly negative 
A for [/ > 2.4. Thus, this phase can be identified with 
the one previously studied^ at A = for n ~ 2.1, which 
is continuously connected with the OSMP at A = and 
half filling {n = 2)^ 

At stronger coupling U > 2.8, additional plateaus ap- 
pear at A < —1.0. In this case, the wide band is half-filled 
and (by the above arguments) insulating. Evidently, the 
crystal field splitting must be essential for this wide-band 
OSMP since it cannot occur at A = (when all phases 
with insulating bands are connected to half filling, i.e., 
particle-hole symmetry). It will, thus, be instructive to 
compare the two different types of OSMPs in some detail, 
in particular regarding spectral properties. 

Before doing that, let us first discuss the impact of the 
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FIG. 2: (Color online) Spectral functions A{u)) for n = 2.1 and 
[/ = 3.0, with data for the wide (narrow) band shown in the 
upper (lower) panels. The QMC estimates (with maximum 
entropy analytic continuation), shown in the left column, in- 
dicate orbital-selective Mott phases of the narrow and (less 
clearly) of the wide band; corresponding Hubbard-I spectra 
capture mainly the peak positions (see text). 



doping level on the basis of Fig.[IJb) which shows results 
for twice the doping strength (filling n — 2.2) compared 
to Fig. [IJa). Overall, the results in both panels of Fig. [1] 
look very similar, up to slight shifts in critical values of 
A and, of course, in the noninteger plateau values. We 
conclude that the specific doping level is not crucial so 
that we can focus on the specific filling n = 2.1 for the 
remainder of the paper without loss of generality. 



tral weight shifts upwards. 

The narrow-band insulating phase is clearly apparent 
as a gap around the Fermi energy for 0.0 < A < 1.2 
in Fig. EJb), in line with the expectations from the or- 
bital occupancy analysis of Fig. [Ija). This gap (with 
a maximum width of about 1.5) shifts with A roughly 
like the center of mass. More generally, the narrow-band 
spectrum starts out from a narrow shape, with a peak 
just below the Fermi energy at A = —3.0, which moves 
towards larger oj with increasing A (and decays slowly 
for A — >■ 3.0). At A « —1.8, a second peak emerges at 
uj w —3, the position of which is initially nearly frozen, 
then moves towards the gap edge at A > 0, until it be- 
comes the main peak. At A = 0.6, the narrow-band spec- 
trum is remarkably symmetric. Minor structures visible 
in the numerical results, such as a splitting of peaks for 
specific values of A, are probably not significant. 

Due to such numerical noise, the wide-band insulat- 
ing phase is much harder to detect in the corresponding 
spectra, Fig.^a). However, a well-developed dip is seen 
to cross the Fermi energy in the range — 1.2 < A < —0.6, 
which is consistent with the very small gap at A —1.0 
that one would expect from interpolating the occupancy 
data of Fig. [TJa). The dip persists at a; « —0.2 for 
A < 1.2, which might be interpreted as a sign of non- 
Fermi- liquid behavior (cf. Ref. |27|). 

In Fig. [2ljc) and Fig. [2jd) we show, for comparison, 
spectra obtained from a Hubbard-I-type solution of the 
same Hamiltonian. Obviously, this simple approxima- 
tion, which docs not include life time effects, cannot cap- 
ture the highly nontrivial correlation physics and repro- 
duce the corresponding spectra on a quantitative level. 
In particular, the Hubbard-I predictions contain unphys- 
ical or much too broad gaps. However, at strong crys- 
tal field splitting |A| — 3.0, the Hubbard-I line shapes 
are roughly correct. In addition, the peak positions and 
their nonuniform evolution as a function of A agree with 
the QMC data remarkably well at all values of A, which 
clearly supports the reliability of our numerical methods. 
For the narrow band, the insulating phase (gap around 
Fermi energy) is also predicted nearly correctly (with a 
false positive only for A — 1.8); in contrast, the wide- 
band insulating phase is grossly overestimated. 



B. Spectral functions 

Now we turn to the spectral properties, starting with 
numerical estimates, obtained by maximum entropy an- 
alytic continuation of QMC imaginary-time Green func- 
tions. Results for the wide band are shown in Fig. HJa) 
and those for the narrow band in Fig. [2fb); correspond- 
ing analytic results (right hand panels of Fig. [2]) will be 
discussed later. As expected, the spectral weight of the 
wide band shifts towards smaller w when the crystal field 
is increased from A = —3.0 (topmost curve in Fig. [2{a)) 
to A = 3.0 (lowest curve in Fig. (^fa)), with a shape that 
changes significantly at intermediate values of A (shown 
with a spacing of 0.6). Conversely, the narrow-band spec- 



C. Phase diagram 

Using primarily orbital occupation data for a broad 
range of parameters A and U , we have constructed the 
phase diagram Fig. [31 Specifically, the critical interac- 
tions for the onset of plateaus, i.e., phases where one of 
the orbitals is incompressible, were determined on a fine 
grid of A values. The resulting boundaries towards these 
orbital-selective insulating phases are shown as symbols 
and solid lines in Fig. [3] for a total filling n = 2.1 (while 
dashed lines denote transitions at n = 2.2). At small 
U ^ 2, the system is either fully metallic or band insu- 
lating (in the narrow band for A —3, in the wide band 
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FIG. 3: (Color online) Phase-diagram as a function of the 
interaction U and the crystal field splitting A for n = 2.1 
(symbols and solid lines) and n = 2.2 (dashed lines). For 
large A, the energetically lower of both orbitals is in a band 
insulating state (BI), while the other orbital contains the re- 
maining electrons. Orbital-selective Mott phases (OSMPs) 
occur at sufficiently strong coupling {U ^ 2) for moderate A. 



for A > 3); the precise critical values of A shift signifi- 
cantly with J7. At J7 > 2, a narrow-band orbital-selective 
Mott phase emerges at A > 0. Only at [/ > 2.8, an ad- 
ditional wide-band OSMP becomes stable at A < —1. 
It is clear that the two OSMPs must be separated (at 
constant filling n > 2) by a fully metallic region, where 
both orbitals are slightly more than half filled [cf. Fig. 
[IJa)]; the extent of this "finger" of the all- metallic phase 
at —1.0 ^ A < —0.2 is practically independent oi U, at 
least in the interaction range covered in this study. 

Stronger doping (n = 2.2, dashed lines in Fig. [3]) sup- 
presses the narrow-band OSMP at A « 0, i.e., shifts the 
phase boundary to larger values of U. Otherwise, the 
effects are surprisingly small, which justifies, again, our 
focus on the single doping level n = 2.1. 



D. Non- Fermi-liquid properties 

In order to get more insight into the nature of the 
orbital-selective phases shown in Fig. |3l beyond the 
orbital-selective incompressibility apparent in Fig. [T] and 
the spectral information of Fig. [51 let us now discuss 
signatures in the self-energy a direct quantita- 

tive measure of correlation effects. Specifically, we con- 
sider the self-energy on the imaginary axis (i.e., for 
u = iujn = i(2n + 1)ttT), which is directly available 
as a state variable in the DMFT self-consistency cycle, 
and avoid an ill-conditioned analytic continuation to the 
real axis (which limits the reliability of spectral functions 
A{uj) = -ImG(w -I- iO+)/7r shown in Fig. [2]). 

The upper panels of Fig. 0] show ImE(ia;„) for the 
wide band (top left) and the narrow band (top right) at 
U = 3.0. Deep in the fully metallic phase, at A = —2.4 
and at A = 3.0 (triangles), both observables have small 
absolute values and decay essentially to zero at small a;„. 
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FIG. 4: (Color online) Imaginary part of self-energy for total 
filling n = 2.1, interaction strengths U = 3.0 (top), U — 
3.6 (bottom), and selected values of the crystal field splitting 
parameter A. 



indicating a weakly correlated Fermi liquid. 

In the narrow-band insulating phase, at A = 0.8 
(squares), the narrow-band self-energy diverges at small 
frequencies; the corresponding wide-band self-energy is 
nearly flat and tends to a finite value at small uj. This 
clear sign of non-Fcrmi liquid behavior in the metal- 
lic component of an orbital-selective Mott phase ap- 
pears completely analogous to the well-known situation 
at n = 2 and A — 0^ Quite remarkably, almost identical 
values of Iml]w(*w„) are obtained also at A = —0.2, very 
close to the edge of the narrow-band OSMP, where the 
narrow-band self-energy ImEn(iw„) is still strongly en- 
hanced at intermediate a;„ ~ 1, but decays almost Fermi 
liquid like at small ujn- 

In the wide-band insulating phase, at A = —1.0 (cir- 
cles), the role of the two bands is just exchanged, relative 
to the case discussed above: now the wide band shows 
a divergent self-energy, while the narrow band displays 
non-Fermi liquid behavior. We conclude that this be- 
havior is really generic of itinerant electrons coupled via 
an Ising type Hund rule interaction to localized electrons 
and does not depend on details of the model. 

A further indication that the behavior discussed so 
far is quite generic is the complete qualitative agree- 
ment between the self-energy data discussed so far for 
U = 3.0 (top panels in Fig. [4]) and corresponding data 
at stronger interaction U — 3.6 (bottom panels in Fig. 
S]): up to a slight enhancement (note the change of scales 
between the two rows of panels), the imaginary parts of 
the self-energies on the Matsubara axis are nearly iden- 
tical at corresponding phase points. The main difference 
is a much larger residual value at small frequencies in 
the fully metallic phases (triangles); apparently, at this 
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stronger interaction (and close to the wide-band OSMP), 
the temperature T = 1/40 is already above the range of 
"good" Fermi liquid behavior at A = —2.4. 



IV. CONCLUSION 

In summary, we have investigated the effects of crystal 
field splitting A in the doped two-band model ([2]) with 
a band width ratio of ^w/^n = 2 and Ising-type exchange 
interaction. Using the Hirsch-Fye quantum Monte Carlo 
technique, we have calculated orbital filling factors, spec- 
tra, and Matsubara self-energies within DMFT, focussing 
on Mott physics within the spatially homogeneous phase. 
The resulting phase diagram contains not only an orbital- 
selective Mott phase of the narrow, i.e., generically more 
strongly correlated, band, but also a wide-band OSMP 
at suitable values of the crystal field (and strong enough 
interaction). This shows that the crystal field, i.e., a 
diagonal element in the tight-binding Hamiltonian, can 
be as relevant for Mott physics as the hopping, i.e., off- 
diagonal elements. Clear signatures of non-Fermi-liquid 
behavior are seen in the Matsubara self-energies in both 
types of OSMPs; these findings are also consistent with 
the complex evolution of the local spectral functions as 
a function of A that we observed using QMC and the 
maximum entropy method. 

On the one hand, the QMC results presented in this 
paper complete the picture regarding the impact of a 
band width difference on the correlation physics of mul- 
tiorbital systems. They show that the OSMT scenario 
established in earlier studies is not only stable with re- 
spect to doping,-- but also with respect to (additional) 
crystal field splitting. On the other hand, our finding of 
a wide-band OSMP makes contact with earlier studies 
of the isolated impact of crystal field splitting on oth- 
erwise equivalent orbitals^"— and shows that the com- 
bined effect of filling control and crystal field^^ can drive 
an OSMT even "against" a significant bandwidth differ- 
ence. 

Our numerical results in the doped case show that the 
range of crystal field splitting A over which the orbital- 
selective Mott phases extend increases significantly with 
interaction U (and Hund rule couplings V = U/2 and 
Jz = U/4). This can also be expected at half fill- 
ing, in line with our analytical finding that the effec- 
tive Heisenberg Hamiltonian relevant in this case and at 
strong coupling does not explicitly depend on A, which 
is, consequently, also true for magnetic ordering tem- 
peratures. It would clearly be interesting to explore 
magnetic order also numerically away from half filling, 
possibly with orbital-dependent frustration.'^^ However, 
quantitative accuracy would then require treatments be- 
yond DMFT)^ at high numerical cost and probably 
with severe sign problems. In contrast, DMFT has been 
shown to be reliable for Mott physics as explored in this 
study; thus, our results are expected to be accurate for 
three-dimensional systems and experimentally relevant. 



e.g., at the lowest temperatures accessible with ultracold 
fermions on optical lattices. 
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Appendix A: Strong-coupling spectral functions at 
half filling 

In this appendix we summarize the calculation of the 
spectral functions in the ground state (T = 0) of the 
two-band model at strong coupling and half filling. The 
derivation proceeds along the lines of Kalinowski and 
Gebhard,'^^ in which the sm(?Ze-band Hubbard model was 
studied in strong-coupling perturbation theory. Here we 
focus on the generic case Ji < Jz, which is somewhat 
simpler and, moreover, is the relevant model for the nu- 
merical calculations in this paper. Calculations for the 
special case J± — Jz — J proceed similarly, in principle, 
but they are technically more involved in detail. 

We start with the definition of the local Green function 
Gmait), which contains (and is in fact equivalent to) the 
local spectral function: 

Here we introduced Heisenberg operators Cj^^(t) = 

e'^'cj^^e"'^* on site 1 for orbital m and spin cr, where 
K ^ H — iiN is the grand canonical Hamiltonian. More- 
over, A/" is the number of lattice sites, T the time ordering 
operator and (• • •)gs an average over all possible degen- 
erate U = oo ground states. The ground state energy at 
half filling is denoted by Eq . The simplest way for impos- 
ing the restriction to half filling in the non-particle-hole 
symmetric two-band model 77a with A 7^ is to con- 
struct a four-band model by duplicating the two-band 
Hamiltonian: 

where the (m = 1,2)- and (m — 3, 4)-orbitals have crys- 
tal field parameters A and —A, respectively. The re- 
sulting four-band Hamiltonian is particle-hole symmetric 
under the transformation — >• (— l)*c|^_|_2 so that 
the chemical potential of the four-band model is exactly 
given by ^ = + U' — ^Jz- Results for the original 
two-band model are then obtained by simply restricting 
consideration to the (m = 1, 2)-orbitals in the end. 

A Fourier transform from the time to a frequency vari- 
able yields two contributions, one from negative and one 



8 



from positive frequencies, which are associated with the 
lower and upper Hubbard bands, respectively: 

The contributions from the Hubbard bands are 

1 



,-n+ '^lm<T)GS 



1 



U ^^''1'"'^ w + K-Eq- iO+ 'i^^-^GS 
and are connected to the spectral function according to 



with ly^y^^iuj) = ilni[G™(w)] for the lower and 
i/UHB(^) ^ -ilm[GUHB(a;)] for the upper Hubbard 
band. It is sufficient to calculate the four contributions 



,LHB 



(cj) to the lower Hubbard band, since the results for 
the upper band then follow automatically from particle- 



hole symmetry, i.e., = iy^^J^ 

In order to calculate we perform a canonical 



transformation to new particles: cj^^ = e^^'^'^c\, 



-S{c) 

which (by definition) leaves the total number of double 
occupancies Hq{c) of these new particles invariant. The 
result is: 



1 



where GS denotes the ground state in terms of the new 
particles, and moreover 



K 



H^ic) + Ho{c) - fiN (Al) 



(ij) 

with m = 1,2,3,4 corresponding to m = 2,1,4,3, re- 
spectively. As a result. 



(ij) 



where the latter term is negligible in the strong-coupling 
limit U/t oo. With the definition 

one then finds that, to this order in the t/[/-expansion. 



simply describes the hopping of a single hole through 
a random environment of dl^\0) and dj^lO) sites. The 
dependence of the right hand side on z already shows 
that the lower Hubbard bands are centered around uj = 
~^{U + Jz ± A), with the sign of ±A depending on the 
band index m. 

In order to calculate G^™(a;), we consider the series 
expansion 



^ oo 



fc=0 



2z[l- S'(z)] ' 



where S{z) describes the contributions to the series of 
hole motions such that the hole does not return to site 
1 between start and finish. The factor ^ in the second 
step occurs, because site 1 has df^jO) -character only with 
probability ^- Furthermore, in each hop only half the 

sites are accessible, because the hole must hop to a dj^|0)- 
site. Hence, on a Bethe lattice (with coordination Z and 



S{z) 



and N = e~^'''^'> Ne^'-''\ Note that the right hand side 
of Eq. (|A1|) defines the effective hopping H^{c) of the 
new particles. We now consider a single hole, i.e., the 
removal of a single particle with spin a from band m, 
in an otherwise half-filled assembly of states dj^jO) and 

dj^jO) with dicr — Ci2crCiicr, which (as argued in SecHI]) are 
lowest in energy ii < J± < Jz- In the following we also yielding 
need single occupancies Sj^^jO) with Sima = Cima- The 
effective hopping H[{c) for a single hole then follows from 
standard Harris-Lange degenerate perturbation theory as 



fc'=0 



[S[z)] 



2z2[l_5'(z)] 



where the sign of the square root is chosen such that 
G™(w) \ for cj -oo. For the (m = 1)- and 
(to = 2)-orbitals, this immediately yields the results for 
the spectral functions quoted in Sec.|lI]for J± < Jz- 

As mentioned above, calculations for the special case 
J± — Jz — J are analogous but considerably more in- 
volved, so that we prefer to simply quote the result (see 
SecHI. 
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